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The Be/X-ray binary A 0535+26 showed a normal (type I) outburst in August/September 2005, 
which reached a maximum X-ray flux of 400mCrab in the 5-lOOkeV range. The outburst was 
observed by INTEGRAL and RXTE. 

The energy of the fundamental cyclotron line has been measured with INTEGRAL and RXTE at 
~45keV. Flaring activity was observed during the rise to the peak of the outburst. RXTE obser- 
vations during one of these flares found the energy of the fundamental cyclotron line shifted to a 
significantly higher position than during the rest of the outburst, where it remains constant. Also, 
the energy-dependent pulse profiles during the flare differ significantly from the rest of the out- 
burst. These differences have been interpreted with the presence of magnetospheric instabilities 
at the onset of the accretion. 

A decomposition method is applied to A 0535+26 pulse profiles. Basic assumptions of the method 
are that the asymmetry observed in the pulse profiles is caused by non-antipodal magnetic poles, 
and that the emission regions have axisymmetric beam patterns. Using pulse profiles obtained 
from RXTE observations, the contribution of the two emission regions has been disentangled. 
Constraints on geometry of the pulsar and a possible solution of the beam pattern are given. First 
results of the comparison of the reconstructed beam pattern with a geometrical model that includes 
relativistic light deflection are presented. 
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1. Introduction 

The Be/X-ray binary A 0535+26 was discovered in 1975 by [[T]]. It is a transient source, in 
an eccentric orbit (e=0.47) of ~ 110.3 days orbital period [^]. Since its discovery it has shown 5 
giant (type II) outbursts in 1980 ig], June 1982 [§, March 1989 [f], February 1994 [g and May 
2005 The last one in 2005 could not be observed in X-rays by most astronomical instruments 
due to Sun constraints. A subsequent normal (type I) outburst took place in August/September 
2005 1^] associated with the periastron passage. The outburst led to our INTEGRAL and RXTE 
observations. 

The source spectrum presents cyclotron resonance scattering features at ~45 keV and ~100 keV 
discovered by HEXE The INTEGRAL and RXTE observations have provided the first high 
accuracy measurements of the cyclotron line energy of A 0535+26. RXTE monitoring has allowed 
to study the evolution of the cyclotron line energy at different luminosities, and like in other sources 



(10], [10], [11]) it was found to vary with the X-ray luminosity [12]. 

Flaring activity has been observed during the rise to the peak of the outburst. RXTE observa- 
tions from one of these flares have revealed a different spectral and timing behavior of the source. 



These results have been presented in [12], and a possible interpretation of the observed changes 
is proposed in []T3[]. A review of the observations and interpretation is given in Sec. ^. In Sec. || 
results from a decomposition analysis of A 0535+26 energy-dependent pulse profiles are shown. 



The method and assumptions are described in Sec.|3. \\ A more detailed description is given in 
Jli]], and therefore all the formal derivations and technical details are not repeated here. We find a 
possible decomposition of the pulse profiles, that has allowed to reconstruct the visible part of the 



beam pattern of A 0535+26, and to extract information on the geometry of the system (Sec. 3.2). 



The reconstructed beam pattern is used as input for a geometrical model that takes into account 



general relativity effects. Preliminary results are presented in Sec. |3.2.3| . Sec. |4] contains a summary 
and conclusions. 



2. IMagnetospheric instabilities 

Flaring activity has been detected during the rise to the peak of the 2005 outburst. RXTE ASM 
and Swift BAT light curves are shown in Fig. |l|. INTEGRAL performed a ~ 200ks observation 
near the peak of the outburst, and RXTE monitored the complete outburst with an exposure time 
of ~ 140 ks, including three observations of one flare during the rise to the peak. The INTEGRAL 
and RXTE observations are indicated in Fig. |I[ 

The pulse period of the pulsar has been studied in detail during the outburst. A constant 
spin period (within uncertainties) is measured during first part of the outburst, P=103. 3960(5) s, 
followed by a spin-up of /"=(— 1.69 ±0.04) x lO^^ss^^ (measured at MJD 53618) at periastron. 
For the first time during a normal outburst a spin-up is measured for A 0535+26, providing evidence 



for the presence of an accretion disk [12]. 



Energy-dependent pulse profiles during the pre-outburst flare differ considerably from those 
obtained during the main part of the outburst, as shown in Fig. g (left). During the main part of the 
outburst, a strong change of the pulse profiles above the cyclotron energy takes place. However, 
during the flare there is a smooth evolution of the pulse profiles above the cyclotron energy. The 
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Figure 1: A 0535+26 RXTE ASM and Swift BAT light curves during the normal (type I) outburst in 
August/September 2005. The dashed vertical line indicates the periastron time. The RXTE and INTEGRAL 
observations are indicated in the upper panel. 



drastic change above the cyclotron energy during the main peak has been interpreted in [ [l3| ] as 
a weakening of the fan-like component of the beam due to the energy dependence of the cross- 
section for ordinary photons. The cyclotron line energy remains constant during the main part of 
the outburst. Near the peak, INTEGRAL measured £'cyc=45. 9 ±0.3keV [15]. A significant change 
in the cyclotron line is measured during the pre-outburst flare, reaching fcyc^Sl.O^j^keV. The 
cyclotron line evolution during the outburst is shown in Fig. ^ (right). 

The flaring activity observed has been explained as due to low-mode magnetospheric instabil- 
ity that develops at the onset of the accretion, in the thin boundary layer between the accretion disk 
and neutron star magnetosphere. The matter accumulated in the boundary layer rapidly falls onto 
the neutron star surface close to the magnetic poles, along different magnetic field lines by which 
stationary accretion proceeds. This explains the differences in the pulse profiles and the change in 
the cyclotron energy during the flare compared to the rest of the outburst. Details are given in [p^]. 
The present analysis will be performed (and is applicable) only to the stable pulse profiles in the 
main outburst. 



3. Pulse profile decomposition 

Geometrical models of filled and hollow accretion columns of accreting neutron stars, includ- 
ing relativistic light deflection, have been computed in [16] and [17]. These models give the beam 
pattern or energy-dependent flux of one emission region as a function of the angle, as seen by a 
distant observer. Introducing the rotation of the pulsar and its geometry, i.e., the orientation of 
the rotation axis with respect to the direction of observation and the location of the two poles, the 
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Figure 2: Left: A 0535+26 energy-dependent pulse profiles during the flare and near the peak of the Au- 
gust/September 2005 normal (type I) outburst. Right: PCA light curve of the outburst (top) and evolution of 
the cyclotron line energy (bottom). 



pulsed emission from each of the two poles or single-pole pulse profiles that a distant observer 
would see can be modeled. The sum of the single-pole contributions gives the total pulse profile. 

In this work we apply a decomposition method that inverts the process described above. From 
the observed pulse profile the single-pole pulse profiles are obtained, and from the single-pole 
contributions information on the geometry of the neutron star and the beam pattern is extracted. 
This method has been applied to the accreting pulsars Cen X-3 [ [T^ ] and Her X-1 [19]. Of course 
the backwards process is not straightforward and involves ambiguities which will be discussed. 



3.1 Method 

3.1.1 Assumptions 

The basic assumption, which is often adopted in model calculations, is that the emission re- 
gions at the magnetic poles are axisymmetric. The beam pattern from one pole is then only a 
function of the angle 6 between the direction of observation and the magnetic axis. This makes 
the single-pole pulse profile necessarily symmetric. One of the symmetry points will be the instant 
when the magnetic axis is closest to the line of sight, and the other symmetry point will be half a 
period later, with the magnetic axis pointing away from the observer. 

Assuming an ideal dipole magnetic field, the above assumption necessarily implies a sym- 
metric total pulse profile. It can be shown that a small displacement of one of the magnetic poles 
from the antipodal position can explain the asymmetry in the sum of the single-pole contributions. 
Therefore the assumption of an ideal dipole field is modified, introducing a small offset from an 
ideal dipole field. 

Another assumption made is that the two emission regions are the same, i.e., have the same 
beam pattern. This implies that each of the two poles will make visible one section of the same 
beam pattern. Depending on the geometry of the neutron star and the angle of observation, those 
two sections will in some cases have coincident parts. This assumption has been tested with the 
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accreting pulsars Cen X-3 [18] and Her X-1 [19|. In those cases an overlapping region was found, 



in agreement with the assumption of two equal emission regions. 

3.1.2 Decomposition into single-pole pulse profiles 

The first step of the analysis is to express the original pulse profile as a Fourier series. The 
total pulse profile F is written as: 

1 "'^-^ n 
F{^) = -uo+ [ukCos{k^)+Vksm{k^)] + u„/2Cos{-^) (3.1) 



where n is the number of bins of the original pulse profile and <I> is the phase. Eq. |3JJ gives 
a valid representation of the original pulse profile at all phases if the Fourier transform of F ap- 
proaches to zero as the frequency approaches n/2. If this is not the case (phenomenon called 
aliasing), less Fourier coefficients are taken into account to describe the original pulse profile. 

The single-pole pulse profiles /i and are described by the following symmetric functions: 

1 n/2 

/i (<!>) = -CO + £c, COS [^(O-Oi)] (3.2) 

1 n/2 

/2(0) = -Jo + £ dkCos{k[^ - (O2 + 71)]} (3.3) 

^ k=l 

<I>i and <I>2 are the symmetry points of fi{^) and /2(<I>) respectively. Formally, a decompo- 
sition of F into two symmetric functions exists for every choice of their symmetry points <I>i and 
<I>2. For convenience, we use the parameter A := tt — (<I>i — O2), that represents the azimuthal 



displacement of one pole with respect to the antipodal position (see Sec. |3.1.3| ). All formal decom- 
positions will be contained in the parameter space — A, with < < 71 and < A < 7r/2. Once 
the formal decompositions are calculated, physical criteria are applied to decrease the number of 
decompositions to physically meaningful ones. Those criteria are the following: 

• "positive criterion": both /i(<I>) and /2(<I>) must be positive, since they represent photon 
fluxes. 

• "non-ripples criterion": the single-pole contributions /i(<I>) and /2(<I') should not be much 
more complicated than the original pulse profile. Individual pulse profiles with many peaks 
that cancel out in the sum are not accepted. 

• the same symmetry points must give valid decompositions in all energy bands. 

Once a possible decomposition is found, the symmetry points for each of the two poles <I>i and 
<I>2, and the parameter A, related to the position of the emission regions on the neutron star, are 
determined. 
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rotation axis 




Figure 3: Schematic view of a rotating neutron star. A spherical coordinate system is used, with the rotation 
axis as polar axis. 0, is the polar angle of the ith pole. 0o is the polar angle of the direction of observation. 
The angle Q between the magnetic pole and the direction of observation changes with the rotation angle <t>. 
Figure from [ |l4] ] . 



3.1.3 From single-pole pulse profiles to beam pattern 

In Fig. ^ a schematic view of a rotating neutron star is shown. A spherical coordinate system 
is used with the rotation axis as polar axis. As was explained above, the beam pattern is assumed 
to be axisymmetric, and therefore to depend only on the angle between the direction of observation 
and the magnetic axis d. The value of 6 changes with the rotation angle 0. Depending on the 
position of the poles with respect to the rotation axis and depending on the direction of observation 
with respect to the magnetic axis, we will only observe a section of the beam pattern for each pole. 

Applying the cosine formula to the spherical triangle in Fig. ^, we obtain as a function of 
the phase 0: 

COS0 =cos0ocos0, + sin0osin0,cos(<l> — <!>;) (3.4) 

where ©o is the polar angle of the direction of observation, 0; the polar angle of the ith pole, 
and <I>, one symmetry point for the ith pole. 

The intrinsic pulsar geometry is shown in Fig. 0. A complete description of the pulsar can 
be given in terms of the polar angles ©i and 02, and the difference in their azimuthal angles 
<I>i — <I>2 = TT — A. The angular distance 5 between the location of the second magnetic pole and 
the point that is antipodal to the first magnetic pole can be used as a measure for the deviation from 
an ideal dipole field. From Fig. ^ 



cos 8 = — cos 02 cos 01 + sin 02 sin 0i cos A (3.5) 

Considering the beam pattern as a function of cos 6 and the single pole pulse profiles as func- 
tions of cos(<I> — <I>,), as the relation between cos0 and cos(<I> — <!>;) is linear (Eq.3.4), there is 
no distortion between the two functions. Therefore, once we have the single-pole pulse profiles, 
plotting them as a function of cos(<I> — O,) we will obtain two sections of the beam pattern. In 
some cases an overlapping region will emerge. Using geometrical properties, it can be shown that 
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2nd Magnetic 
Pole 



Figure 4: Intrinsic geometry of the pulsar. With the rotation axis as polar axis, the magnetic poles are 
located at polar angles ©i and 02. The angular distance 5 between the second magnetic pole and the point 
that is antipodal to the first magnetic pole gives the deviation from an ideal dipole field. Figure from [p^. 

the relative position of the two single-pole profiles (in terms of the total pulse profile) are related 
to the position of the two poles. This requires an independent determination of the direction of 
observation. Details can be found in [pl|]. 

3.2 Application to A 0535+26 

RXTE HEXTE energy-dependent pulse profiles of A 0535-1-26 during the August/September 
2005 normal outburst have been analyzed. Pulse profiles obtained during the main part of the 
outburst have been selected for the analysis, as those profiles appear to be very stable, not only 
during the outburst but also compared to historical observations. Several observations have been 
used, for a total exposure time of ~20ks. 

3.2.1 Search for acceptable decompositions 

The energy-dependent pulse profiles are written as Fourier series. These functions are then 
written as the sum of two symmetric functions /i(<I>) and /2(<I>). To search for physically mean- 
ingful decompositions, the — A parameter space is divided in 1 ° x 1 ° boxes. All the formal 
decompositions are represented in this plane. We apply the "non-negative criterion". Out of all the 
decompositions, only those for which /i (O) and /2 (O) are positive are accepted. Fig. ^ (left) shows 
where in the — A plane positive decompositions have been found. The "non-ripples criterion" 
is then applied. A quality function is defined, which counts the number of peaks in one profile. To 
avoid having to study all the decompositions, similar ones are grouped into types. We study the 
representative of each type. Fig. |5| (right) shows where the highest ranked profile representatives 
were found. We have examined all these decompositions in different energy ranges, dividing the 
parameter space in five regions A, B, C, D and E. 

Combining the results of the analysis in different energy ranges and for different observations, 
we find a best decomposition of the original pulse profiles in region A. They are shown in Fig. ^. 
The dotted horizontal line corresponds to the remaining unmodulated flux available to distribute 
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Figure 5: Left: Result of applying the "non-negative" criterion: decompositions outside the black area are 
discarded. Right: Highest ranked decompositions grouped in five regions A, B, C, D and E. The different 
symbols represent different energy ranges: +: (18.3-30.9)keV, x: (30.9^4.5)keV , A: (44.5-59. l)keV, 
□: (59.1-99.8)keV 
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Figure 6: Best decomposition of the original pulse profiles (left) in two symmetric functions (right). 



between the two symmetric functions that can not be determined from the decomposition. The 
minima of the two symmetric functions have been shifted to zero, so that the sum of the two 
symmetric functions plus the unmodulated flux reproduces the original pulse profile. The values 
of , A and <I>2 (using A = k — (Oi — O2)) can not be more accurate than the 12° bin resolution, 
and we therefore use this value as uncertainty for the symmetry points. Average values for the best 
decomposition in region A are = 73° ± 12°, A = 33° ± 12° and + TT = 106° ± 12° 

Decompositions in regions C and E are discarded because they present a very strong anti- 
correlation in the main peaks which seems artificial, not expected from two independent emission 
regions. Decompositions in regions B and D have also been discarded, because the single-pole 
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Figure 7: Reconstructed beam pattern of A 0535+26 in different energy ranges. Left: linear representation. 
Right: polar diagram. 



pulse profiles present an anti-correlation in many small features that cancel out in the sum, also not 
expected from two independent emission regions. Another argument to reject decompositions in 
regions C, D and E is that they all present higher values of A. Under the assumption of slightly 
displaced magnetic poles, smaller values of A are more likely to be real. This was the case in the 



analysis of the accreting pulsars Cen X-3 [18] and Her X-1 [19], where the best decompositions 
were found for small values of A. A further argument against decompositions in regions B, C 
and D emerges in the reconstruction of the beam pattern from the single-pole contributions (see 



Sec. 3.2.2). 



3.2.2 From single-pole pulse profiles to geometry and beam pattern 

Doing the proper transformation of the single pole pulse profiles into undistorted sections of 
the beam pattern, an overlapping region does not emerge. However, under the assumption of two 
identical emission regions, the two sections can almost be connected to each other, with a small gap 
in between. This was not possible for the decompositions in regions B, C and D. The shift between 
the two sections is related to the position of the magnetic poles and the direction of observation 



(see [ ]14[ ] for details). The inclination of the system is / = 37 ±2° [|20|]. Assuming that the rotation 
axis of the neutron star is perpendicular to the orbital plane, / = 0o- We can therefore obtain the 
location of the poles ©i and ©2- The angular distance between the location of the second pole and 



the point that is antipodal to the first pole 5 can be estimated using Eq. 3.5. The estimated values 
for the position of the magnetic poles and the offset are 0i k, 50°, ©2 ~ 130°, 5 25°. It is then 
possible to plot the reconstructed sections of the beam pattern as a function of 6. The two sections 
of the beam pattern are reconstructed for Q € (13° — 87°) and Q € (93° — 167°). The reconstructed 
beam pattern is shown in Fig. 0. 

3.2.3 Interpretation of beam pattern 

A characteristic feature of the reconstructed beam pattern is a minimum observed in the flux 
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between d ^ 30° —40°. This feature is present at all energies and we therefore believe that it could 
be related to the geometry of the accretion. (It is also present at lower energies, seen from the 
analysis of PCA data [23]). Filled column models give a beam pattern in which the flux decreases 
at low values of 6, corresponding to the instant when the observer looks along the accretion stream 
[17]. Introducing a hollow column plus a halo created on the neutron star surface around the 
column walls from scattered radiation emitted from the walls would explain the flux at low values 
of 6 and the minimum as d increases when the observer looks directly into the column. 

The steep increase in flux at high values of theta {6 > 120°) could be due to gravitational light 
bending, which produces a similar feature in model calculations. 

To obtain estimates on the size of the accretion column, a model for a hollow column plus 
a halo has been calculated for 6 G (0°,40°). Geometrical models of filled columns, including 
the formation of a halo around the accretion column, were presented in [17], where the relative 
importance of the different components (halo-column) to the observed flux was studied. In the 
work presented here, the modelling is performed as in |[l7|], but introducing a hollow column. A 
detailed study of this model will be presented elsewhere. Beam patterns are computed using ray- 
tracing [21] and include relativistic light deflection [22]. 

The emission of the column wall is assumed to be a blackbody of temperature T^ff, assumed 
to be isotropic. Part of the radiation that leaves the column hits the neutron star surface, creating a 
luminous halo around the column walls. This emission is assumed to be thermal. Schwarzschild 
metric is used to calculate photon paths. This is appropriate for A 0535-1-26 because it is a slowly 
rotating neutron star {Pspin ~ 103.4^). The Schwarzschild radius is = 2GM/c^ and the neutron 
star radius r„. The accretion funnel has an inner half-opening angle a, and outer half-opening 
angle ao- A radiative shock is assumed to form close to the neutron star surface at radial coordinate 
Vf. See Fig. || for a sketch of the geometry of the hollow column. Below the shock, the column 
is optically thick and radiation is emitted from the inner and outer walls. Above the shock, the 
accretion is assumed to be in free-fall. 

The mass and radius of the neutron star are taken to be M„ = 1.4M0 and r„ = 10km. The 
radiative shock is formed at rt = 10.5 km, and an asymptotic luminosity per pole of Leo = 0.8 x 
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Figure 9: Beam pattern models for hollow column with halo for different column thickness and the same 
opening angle (left) and different opening angles with same column thickness (right), calculated for Q G 

(0°,40°). 



lO^^ergs^' is taken. The cyclotron energy, used to calculate the magnetic scattering cross-section, 
is set to £'cyc=45 keV at the magnetic pole. Beam patterns have been computed for 7.6 keV photons. 

Different models have been computed for different values of the inner and outer half-opening 
angles. For each model, the temperature of the Plank spectrum emitted by the column wall and the 
density that the accreting material has at the base of the free-fall section follow from the assumed 
geometry and the value of Leo. The parameters used for each model are listed in Table |l|. Results 
from the model calculations are shown in Fig. ^ 

Table 1: Model parameters used in the computation of beam patterns for A 0535+26, for Q E (0°,40°). 
model a, (rad) ao (rad) ^r(keV) p (lO^^g/cm^^) 



1 


0.08 


0.1 


4.1 


16 


2 


0.06 


0.1 


4.1 


9 


3 


0.04 


0.1 


4.1 


6.8 


4 


0.09 


0.15 


3.7 


4 


5 


0.14 


0.2 


3.5 


2.8 



By comparing the A 0535-1-26 reconstructed beam pattern (Fig. |7|) with the models in Fig. ^, 
it can be seen that the shape of the models reproduce well the shape of the reconstructed beam 
pattern. We can estimate the half-opening angle and column thickness to be ao = 0.2rad~ 11.5° 
and Oq — at = 0.06rad~' 3.4°. Those are preliminary results from an ongoing work that will be 
presented in [^]. 



4. Summary and conclusions 



In this work results from spectral and timing analysis of A 0535-1-26 during a flare that took 
place in the rise to the peak of the August/September 2005 normal outburst compared to the rest 
of the outburst have been summarized. The differences observed in the pulse profiles and in the 
cyclotron energy can be explained with magnetospheric instabilities that develop in the onset of the 
accretion. Those results and interpretation can be found in [12] and 1 13] respectively. 



11 



A0535+26 in outburst 



I. Caballero 



In the second part of the work, a decomposition analysis has been applied to the A 0535+26 
energy-dependent pulse profiles. A dipole magnetic field is assumed with axisymmetric emission 
regions. The asymmetry in the total pulse profiles is explained with a small offset from one of 
the magnetic poles from the antipodal position. We find a physically acceptable decomposition 
of the pulse profiles that allows to extract information on the geometry of the pulsar. We obtain 
©1 ?s 50° and ©2 ~ 130° for the position of the magnetic poles, and an offset of 5 f« 25°. The 
visible section of the beam pattern has been reconstructed. First results from a geometrical model 
of hollow column that includes the formation of a halo are used to interpret the main features of 



the reconstructed beam pattern of A 0535+26. Further details will be presented in [23]. 
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